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We revisit the transverse single spin asymmetry in the angular distribution of a Drell-Yan dilepton 
pair. We study this asymmetry by using twist-3 collinear factorization, and we obtain the same 
result both in covariant gauge and in the light-cone gauge. Moreover, we have checked the electro- 
magnetic gauge invariance of our calculation. Our final expression for the asymmetry differs from 
all the previous results given in the literature. The overall sign of this asymmetry is as important 
as the sign of the Sivers asymmetry in Drell-Yan. 

r^, 1 1 Introduction 

The observation of transverse single spin asymmetries (SSAs) in various hard scattering processes has 
stimulated new remarkable developments both on the theoretical and the experimental side. As a 
f"*-. ! consequence, the study of SSAs currently represents a very active field of research [THS]. The interest 
in such effects is essentially twofold: first, SSAs allow one to address the parton structure of the 
nucleon beyond the collinear parton model approximation. Second, SSAs are ideal observables in 
order to further explore in which cases the machinery of QCD factorization still applies and in which 
cases, in its simplest form, it breaks down (see [I] and references therein). 

For what concerns the parton structure of the nucleon, in the present work we focus on collinear 
twist-3 quark-gluon-quark correlations. To be more precise, the central non-perturbative correlator 
is the so-called ETQS (Efremov-Teryaev-Qiu-Sterman) matrix element [5HZ] Tp — and its chiral-odd 
partner Tp — which typically appears when describing transverse SSAs in the context of collinear 
higher-twist factorization. The machinery of collinear twist-3 factorization was pioneered already in 
the early 1980's 0[8l[9], and in the meantime frequently applied to transverse spin effects in hard 
semi- inclusive reactions (see, e.g., Refs. [6ll7l [l0l - [T2] ) . 

In this paper, we revisit the transverse single spin asymmetry in the angular distribution of a 
Drell-Yan dilepton pair. This asymmetry is defined as the difference of two spin dependent cross 
sections with opposite directions of transverse polarization divided by their sum, 



in 



_ ( da{S T ) da{-S T ) \ /f da(S T ) da(-S T ) \ 
N V dVLdQ 2 dttdQ 2 J I \ dttdQ 2 dQdQ 2 J ' 



where dO, = dcos 8d(j)s is a solid angle element of the leptons in a dilepton rest frame, and the azimuthal 
angle 4>s is measured relative to the transverse spin vector. Note that the transverse momentum Qt of 
the dilepton pair is integrated out, and we emphasize that integrating over Qt is essential for applying 
the collinear factorization approach in the present case. 

The asymmetry An was already studied in several previous articles, and various different results 
were obtained. The first calculation, carried out in the light-cone gauge, can be found in Ref. |13j . 
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The authors obtained^ 

AHTS) 1 sin20sin05 E g e g I dx (Tp(x,x) — x-^Tp(x,x)^ff(x ) 
N = ~Q 1 + cos^ E^/efc/fMAV) ' 

where is the standard unpolarized twist-2 parton distribution for quark flavor q. The momentum 
fraction x' is given by x' = Q 2 /(xS), with S = (P + P) 2 denoting the square of the cm energy of the 
process. Later on the presence of the derivative term in the numerator of fl2J) was doubted, and it was 
argued that the correct result for An should be [T^[T5] 

Abmt) = _ 1_ smjgsin^g Eg e 2 q J dxT F (x,x)ff(x') 

N Q l + cos'fl E q e 2 q fdxf«(x)f?(x>) ' U 

Afterwards, in Ref. [16] ^4 at was considered in the collinear twist-3 approach, and the result of that 
study agreed with the expression in ([3"j). Then An was computed by using factorization in terms 
of transverse momentum dependent correlators [17J. The final outcome of that work neither agreed 
with ([2]) nor with (|3|) . More recently, An was again considered in Ref. [18] , where the authors claimed 
that the spin-dependent hadronic tensor should be multiplied by a factor of 2 compared to previous 
work [13H16] . The controversy about the derivative term was not addressed in [18] , 

This somewhat unclear situation motivated us to revisit this topic. We computed An in ([T]) by 
means of twist-3 collinear factorization and came up with yet another result. To be specific, our result 
is just half of the one quoted in Eq. ([3]). An important difference in comparison to previous work 
is that, when performing the collinear expansion in the twist-3 formalism, we take into account the 
dependence on transverse parton motion (fey-dependence) not only in the hadronic tensor but also in 
the lepton tensor. In order to gain further confidence we checked our calculation in a few different 
ways. 

The rest of the paper is organized as the follows. In the next section, we introduce our notation, 
and give some details about the kinematics. In Section 3, we derive the asymmetry in a covariant 
gauge as well as in the light cone gauge, and the two results agree with each other. In addition, we 
have checked the electromagnetic gauge invariance by explicit calculation. We summarize the paper 
in Section 4. 



2 Kinematics and notation 

We focus on lepton pair production in hadronic scattering which comes from the decay of a virtual 
photon, Ha + Hb — > 7* + A" — > £ + + A. The 4-momenta of the leptons are l\ and I2, and q = h + h 
denotes the momentum of the virtual photon. The invariant mass of the dilepton pair is Q with 
Q 2 = q 2 . For the following calculation we need to introduce the vector R = l\ — £2- In any dilepton 
rest frame, R reads 

R = Q (0, sin 9 cos <p, sin 9 sin <p, cos 9^j , (4) 

where the numerical values of 9 and <p depend on the frame. The correlation associated with An in ([1]) 
is e flupa P fl P u S p R a , while the asymmetry usually associated with the Sivers effect [19] is related to 
the correlation e iiVpa P lL P v S p q a . The latter requires to measure the transverse momentum Qt of the 
dilepton pair. 

lr To shorten the notation we suppress throughout terms where quarks and antiquarks are interchanged. 
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A convenient way of sorting out the different angular dependences of the Drell-Yan cross section 
is to decompose the lepton tensor in terms of individual independent orthogonal tensors |20H22j . 



= ((q + RY(q - Rf + {q + Rf(q - R)" - 2Q 2 g^\ = ]T L { V^ . (5) 



(A discussion of the general structure of the polarized Drell-Yan cross section can be found in Ref. |23j.) 
The Li represent the angular structures, and the basis tensors can be constructed from a set of 
(4-dimensional) basis vectors T^, X^, Y^, Z^, which are mutually orthogonal to each other and are 
normalized according to T 2 = 1, X 2 = Y 2 = Z 2 = — 1. For the case of An, the relevant angular 
structures appear in the terms associated with V 3 = -\{Z*X V + Z V X> X ) and V s = -\{Z^Y U + Z U Y^), 

If = Q 2 sin 29 cos + Q 2 sin 26 sin Vg v + . . . . (6) 

Like in the case of the decomposition Q holds in any dilepton rest frame. 

For Qt = (in the hadronic cm frame), we choose the following dilepton rest frame: z-axis along 
the direction of the polarized hadron, and x-axis along the direction of the polarization vector St- To 
be fully specific, the 4-dimensional basis vectors are given by 

= 



2 



VQ 

Z^ = — (xP^ - x'P>) , 
X» = S T , 



Y » = e » vpCT T v Z p X a . (7) 

Because of the specific definition of Z^ , this frame can actually be considered as partonic cm frame. 
If Qt 7^ 0, one may work in the Collins-Soper frame [24J for which the basis vectors read 



Z» = 1 2 (qpP^-qpP"), 



Y n = £ ^P°T v Z p X a . (8) 

In © we use the further definitions = [P» - (P • q) / q 2 'q») / 'V~S , P* = [P> - (P • q)/q 2 q> l ]/VS, 
with q p = P ■ q/V~S, qp = P ■ q/y/S. At tree level, Qt is equal to the sum of the intrinsic transverse 
momenta of the two incoming partons. This implies that for Qt / Oa /cr-dependence is sitting in the 
unit vectors X 1 - 1 , Y^ and T M . (The fey-dependence of Z^ is of the order k T and therefore irrelevant 
for our twist-3 calculation.) As a result, the terms containing cos</> and sin</> are /c^-dependent. This 
fey-dependence must be taken into account when performing the collinear expansion. 

For Qt ^ 0, instead of using the Collins-Soper frame, one can alternatively perform the calculation, 
for instance, in the Gottfried- Jackson frame |25j . Keeping track of all fey-dependent terms in the 
Gottfried-Jackson frame is more involved. Nevertheless, we carried out the calculation, and our final 
result agrees with what we find in the Collins-Soper frame. 
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3 Calculation in twist-3 collinear factorization 



In order to calculate A jy in Eq. ([T]) one needs both the unpolarized cross section (in the parton model) 
and the spin-dependent cross section. The former is well-known and given by 

d^n = ^ ? e « / dx dx ' fKx) /lV) [fit (1 + cos2 e) 5 ( q2 - xx ' s )} ■ (9) 

The polarized cross section is a twist-3 effect and depends on quark-gluon-quark correlations, which 
contain interesting physics beyond the parton model. In fact, such twist-3 correlations associated with 
both hadrons can give rise to An leading to the generic expression 



. 1 T F (x,x)®f 1 (x , ) + h 1 (x)®Tp*\x',x') 

Anoc q m^w) ' (0) 

where hi is the transversity distribution. The second (chiral-odd) term in the numerator, which we 
have not included in Eqs. ([2]) and ([3]), was first considered in Ref. [14] . In our calculation we treat 
both the chiral-even and the chiral-odd contribution to An- 

The ETQS matrix element Tp and its chiral-odd partner Tp^ are defined asH 

T F (x, Xl ) = J ^^ e -^V-+i(*i-*)^ 

4 a \x, Xl ) = |^^ e -^v-+*>i-*)^ (11) 

where a summation over color is implicit, and gauge links has been suppressed. In the following two 
subsections we compute the hard coefficients associated with these matrix elements both in covariant 
gauge and in the light-cone gauge. It is worthwhile to mention that Tp(x,x) and T F (x,x) are 
related to particular fcy-moments of the transverse momentum dependent Sivers function [19] and 
Boer-Mulders function [26], respectively [TTtl2"T]. 



3.1 Asymmetry derived in covariant gauge 

In covariant gauge, the leading contribution of the gluon field is from the component parallel to 
the direction of its momentum. If one considers P + (with P being the momentum of the polarized 
nucleon) and P~ as the large light-cone momenta, then the dominant component of the gluon field 
for the diagrams shown in Fig. Q]is A + . Before making the collinear expansion, the incoming partons 
carry a transverse momentum kiT, which is much smaller than the dominant longitudinal momentum. 
In order to extract the twist-3 contributions from the diagrams with one-gluon-exchange, one needs to 
get one power of kiT from the hard scattering part and combine kiT with A + in order to convert the 
gluon field in the matrix element into the corresponding part of the field strength tensor (8tA + ) [7]. 
As stated above, the /c-r-flow may go through the lepton lines via the virtual photon. Therefore, we 
have to expand the hadronic tensor as well as the lepton tensor in terms of kiT around k^ = 0. 

Only the two diagrams in Fig. [1] contribute to An in covariant gauge. To be more precise, these 
two diagrams provide the chiral-even Tp part of the asymmetry. In order to get the chiral-odd 
contribution one has to consider the corresponding two diagrams for which the gluon is associated 
with the unpolarized hadron. As an example, for the left cut-diagram in Fig. Q] we have the following 

2 For a generic 4-vector v, we define light-cone coordinates according to v = (v° ± i> 3 )/\/2 and vt = (v 1 ,v 2 ). 
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Figure 1: Diagrams contributing to An in covariant gauge. The gluon attached to the hard scattering 
part is longitudinally polarized. In order to extract the twist-3 contribution, one has to expand in kx 
and kiT, and to pick up the linear terms. 



expansion, 



H^'P(xp + k T , x lP + k 1T , S T ) P p L^q = x x p + k XT + 



x' 



■'p,R) 



H^'P(xp, xip, S T ) P p Lp V {q = x lP + x'p, R) 

' d ( cos <j> Vg$, H^p P p ) d ( sin 4> Vg£ H^'P P f 



+Q z sm26 



+Q 2 sm26 



■ l rrCMdH^'PPp 

sm<p s V 8 CA 



+ 



Ok 



IT 



K IT 



dkS 



krp + ... , 



(12) 



the tensors V^ v and V^ v , and the hadronic 



where the superscripts CM and CS refer to the partonic cm frame and the Collins-Soper frame 
specified in (|7|) and in (JSj) , respectively. The azimuthal angle <f> is understood in the Collins-Soper 
frame, while the azimuthal angle in the cm frame is just what we defined above as 4>s, namely the 
angle between Rt and St- There is no need to distinguish between the polar angle 6 in the two frames 
when expanding around kiT = and keeping only the linear terms. For the left cut-diagram in Fig. [H 
the lepton tensor is independent of kx, but it depends on k\T- The used tensor decomposition of the 
lepton tensor is rather convenient in order to treat this fciT-dependence. This dependence is sitting 
in three parts: the angular dependences cos 4> and sin 
tensor H^'f P p . 

The first term of the Taylor expansion in (I12D corresponds to the eikonal line contribution to the 
twist- 2 quark distribution, which does not contribute to the asymmetry. One can extract the desired 
twist-3 term by picking up the terms linear in kT (and kvr) from the above expansion. Note that in 
HP U,P P p also a delta function of the form 5(Q 2 — (xp + k\T + x'p) 2 ) is hidden. It is easy to see that 
this delta- function cannot provide a term linear in k\x, and therefore its fciT-dependence is irrelevant 
for the calculation of A^. This is actually the reason why the derivative term of Tp, which we briefly 
discussed in the Introduction, does not show up in A^. In general, the collinear expansion enables 
one to integrate out three of the four components of the parton loop momenta, and as a result the 
non-perturbative part can be expressed through the collinear twist-3 correlations Tp and Tp . 

The strong interaction phase necessary for having a nonzero SSA arises from the partonic scattering 
amplitude with an extra gluon. As is evident from the diagrams in Fig. HJ this amplitude interferes 
with the real scattering amplitude without a gluon. The imaginary part is due to the pole of the quark 
(antiquark) propagator and arises when integrating over the longitudinal gluon momentum fraction 
x g . In the present case, one has a pole for x g = ("soft gluon pole" from initial state interaction), 
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xP + k T 




Figure 2: Contribution from the ^-expansion in the light cone gauge. The /c^-flow goes also through 
the lepton lines via the virtual photon propagator. 



while there is no contribution from so-called hard gluon poles or soft fermion poles. We extract the 
imaginary part of the pole by using the formula 



Im ■ 



x g db ie 



Collecting all the pieces we finally arrive at the following polarized differential cross section, 



dQ 2 dn 



9Q 2 



E 



e\ I dxdx' (T F {x,x)ff(x') + h q 1 {x)T ( p )q (x',x' 



1 

Q 



— (- sin 26 sin (f>) 5(q 2 - xx'S 

327T V 



This provides the asymmetry 



A N = - — 



1 sin 2(9 sin 05 £<? 4 I dx (t f (x, x) ff(x') + h\{x) T£> q (x', x'\ 



2Q l + cos 2 # E q e 2 q Jdxf!(x)f!(x') 
which, as already stated above, is just half of the result ([3]) obtained in Refs. [14TH6] . 



(13) 



(14) 



(15) 



3.2 Asymmetry derived in the light-cone gauge 

To test the color gauge invariance of our result, we derived the asymmetry also in the color light-cone 
gauge. In general, in the light-cone gauge both the first order /cy-expansion of the born diagram 
(see Fig. [2]) and the diagrams with one additional exchange of a transversely polarized gluon (see 
Fig. [3]) contribute to the spin dependent cross section at the twist-3 level. The associated twist-3 non- 
perturbative parts are the matrix elements for which the operators tpdTip and i/jAtiP are sandwiched 
between the hadron state [8]. Apparently, these two correlators are not QCD gauge invariant. However, 
if one entirely fixes the light-cone gauge, i.e., if one carries out the calculation using a specific boundary 
condition for the transverse gluon field at the light-cone infinity, then the two matrix elements can be 
uniquely related to the gauge invariant quark-gluon-quark correlators Tp and Tp^ |12|I29]. 

There exist three frequently used boundary conditions: the retarded boundary condition, the 
advanced boundary condition, and the anti-symmetric boundary condition. For the Drell-Yan process 
the retarded boundary condition At{— oo~) = is the most convenient choice [30J. Exploiting this 




(a) 



(6) 




(c) 



Figure 3: Feynman diagrams with one-gluon-exchange relevant for the calculation of A jy in the light- 
cone gauge. The momenta carried by all incoming partons have only a longitudinal component. 
The diagrams (c) and (d) represent the contribution from the so-called special fermion propagator 
introduced in Ref. [28]. Note that the special propagator actually contributes to the hard coefficients 
of both the operator i/jOt^ and the operator ipATip- One finds that these two contributions exactly 
cancel each other. 



particular boundary condition, the operators ipdTip and ipAxip can be readily rewritten in a gauge 
invariant form. For example, one has 



dy 
An 



JxP+y 



(PS\iP(0) 7 + 4" Stu i&r» i/>(y~)\PS) = T F (x, x) 



(16) 



as well as 

dy'dy^ 



Air 



P +^P + y-^-*i)P + yT (PS^(O) 7 + S Tu gA T ^) ij>(y-)\PS) 



dy dy 



An 



1 e ixP+ 



y-^- Xl )p +y - {PS $(0) 7 + £ v / Stu gF + (y^)4,{y-)\PS) . (17) 



X — X\ + l€ 

One has to organize the contributions associated with ipdr^ and i/jAtiIj i n a different way when using 
different boundary conditions [12,29]. Though the final result is of course independent of the boundary 
condition, the calculation of the hard part associated with tpA^ in Drell-Yan is much more involved 
for both the advanced and the anti-symmetric boundary condition. For a general discussion about 
these issues and some more technical details we refer the interested reader to a forthcoming paper [31] . 
For the chiral-even contribution, the generalized factorization formula takes the form 



da(ST) 



dQ 2 dn 



12Q 2 *f q 



(~ I dx dx' Tp{x, x) fi(x') e^f St p 



_d 

dk: 



V ( H Bom(xP + kr, x'p) V 3 ™ sin 29 COS , 
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+ H Zrn( x P + k T, x'P) V 8 C ^ sin 26 sin 
1 f i 

H — / dxi H^(xp,xip,x'p)Vf.^ sin 20 sin </> s 

7r J x — xi+ie 



k T =0 

(18) 



where, in the end, only the /c^-expansion of the hadronic tensor contracted with the tensor V^ v 
contributes to the asymmetry, while the corresponding expression associated with V%^ v vanishes due 
to parity conservation. This point is exactly reversed in the case of the chiral-odd part related with 
Tp\ Note that the required imaginary part in the hard term coupled with the operator •ipATip can 
arise from the (artificial) pole l/(x — xi + ie), which is generated by partial integration in Eq. (|17p. 
Moreover, the diagrams with a special propagator contribute to the hard parts resulting from both 
the /^-expansion and the gluon-exchange. However, these two contributions cancel each other. 

The perturbative calculation is rather straightforward. The final result for An of the calculation 
in the light-cone gauge exactly matches with the final result (|15p we found in covariant gauge. 



4 Summary 

In summary, we recalculated the transverse single spin asymmetry An in the angular distribution 
of a Drell-Yan dilepton pair by using twist-3 collinear factorization. Compared to previous work on 
this topic, we payed particular attention to the fc^-dependence of the lepton tensor when making the 
collinear expansion. Our final result for An in Eq. (I15j) differs from all the previous results given in the 
literature. For instance, we find an asymmetry which is just half of what was obtained in Refs. [14H16] . 

We made various checks in order to gain further confidence in our calculation. First, we verified 
QCD gauge invariance by performing the calculation both in covariant gauge and in the light-cone 
gauge. Second, we tested the electromagnetic gauge invariance by recalculating the asymmetry in 
two specific QED light-cone gauges. Third, we computed the NLO real emission corrections in the 
leading-log approximation. In a certain sense, this calculation is more straightforward than the lowest 
order treatment, since the /c^-flow can go through the unobserved parton line. The outcome of this 
study fully supports our result for An presented in the present work. A complete NLO analysis will 
be presented elsewhere. 

It is important to notice that measuring the sign of An can be considered to be equally important 
as checking the predicted sign reversal of the Sivers effect in Drell-Yan [32] . In either case the physics 
of initial state gluon interactions would be tested. 

The formalism developed in this paper can be extended in order to study similar observables 
which represent a correlation between the transverse spin and the relative transverse momentum of 
final state particles. For instance, transverse SSAs for dihadron production in semi-inclusive DIS can, 
in principle, be treated along the same lines. We plan to address this point in a future work. 

Acknowledgements: We thank F. Yuan and D. Boer for helpful discussion. This work is supported 
by the NSF under Grant No. PHY-0855501. 
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